“Calhoun 


Institutional Archive of the Naval Postgraduate School 





Calhoun: The NPS Institutional Archive 
DSpace Repository 


Theses and Dissertations 1. Thesis and Dissertation Collection, all items 


2001-09 


A nonlinear wave shoaling model for 
alongshore varying Bathymetry 


Ruth, David M. 


Monterey, California. Naval Postgraduate School 
http://hdl.handle.net/10945/1450 


This publication is a work of the U.S. Government as defined in Title 17, United 
States Code, Section 101. Copyright protection is not available for this work in the 
United States. 


Downloaded from NPS Archive: Calhoun 


Calhoun is the Naval Postgraduate School's public access digital repository for 


\§ D U DL EY research materials and institutional publications created by the NPS community. 
«iis Calhoun is named for Professor of Mathematics Guy K. Calhoun, NPS's first 


NNN KNOX appointed -- and published -- scholarly author. 


LIBRARY Dudley Knox Library / Naval Postgraduate School 
411 Dyer Road / 1 University Circle 


http://www.nps.edu/library Monterey, California USA 93943 


NAVAL POSTGRADUATE SCHOOL 
Monterey, California 





THESIS 


A NONLINEAR WAVE SHOALING MODEL FOR 
ALONGSHORE VARYING BATHYMETRY 


by 


David M. Ruth 


September 2001 


Thesis Advisor: Thomas H.C. Herbers 
Second Reader: Edward B. Thornton 





Approved for public release; distribution is unlimited 


Public reporting burden for this collection of information is estimated to average 1 hour per response, including 
the time for reviewing instruction, searching existing data sources, gathering and maintaining the data needed, and 
completing and reviewing the collection of information. Send comments regarding this burden estimate or any 
other aspect of this collection of information, including suggestions for reducing this burden, to Washington 
headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 
1204, Arlington, VA 22202-4302, and to the Office of Management and Budget, Paperwork Reduction Project 
(0704-0188) Washington DC 20503. 


September 2001 Master’s Thesis 
A nonlinear wave shoaling model for alongshore-varying bathymetry 
7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 8. PERFORMING 


Naval Postgraduate School ORGANIZATION REPORT 
Monterey, CA 93943-5000 NUMBER 


9. SPONSORING / MONITORING AGENCY NAME(S) AND ADDRESS(ES) 10. SPONSORING / MONITORING 
N/A AGENCY REPORT NUMBER 


11. SUPPLEMENTARY NOTES The views expressed in this thesis are those of the author and do not reflect the official 
policy or position of the Department of Defense or the U.S. Government. 


12a. DISTRIBUTION / AVAILABILITY STATEMENT 
Approved for public release; distribution is unlimited. 

13. ABSTRACT (maximum 200 words) 

This thesis proposes an improvement to present near-shore wave prediction models. Using 
weakly dispersive Boussinesq theory, the shoaling of directionally spread surface gravity 
waves over a beach with gentle gradients in the cross-shore and alongshore directions is 
examined. Following Herbers and Burton (1997), the governing fluid flow equations are 
expanded to third order and depth-integrated over the water column. A resulting amplitude 
evolution equation for a spectrum of waves is derived, which is the main result of this paper. 
New terms in the higher order result include effects due to alongshore bottom slope, higher 
order cross-shore depth variations, and non-linear quartet interactions. The linear terms in this 
equation are verified by analytical methods using linear finite depth theory. Example 
computations for a monochromatic wave train over a plane beach quantify some of the 
improvements of this result over the lower order model. Opportunities for further development 
and verification of this result are proposed, and recommendations for application of the result 
in its present form are outlined. 


14. SUBJECT TERMS Boussinesq Wave Model, Ocean Surface Gravity Waves, Surf Zone, Wave | 15. NUMBER OF 


Shoaling, Beach. PAGES - 


16. PRICE CODE 


17. SECURITY 18. SECURITY 19. SECURITY 20. LIMITATION 
CLASSIFICATION OF CLASSIFICATION OF THIS CLASSIFICATION OF OF ABSTRACT 
REPORT PAGE ABSTRACT 
Unclassified Unclassified Unclassified UL 
NSN 7540-01-280-5500 Standard Form 298 (Rev. 2-89) 
Prescribed by ANSI Std. 239-18 





THIS PAGE INTENTIONALLY LEFT BLANK 


il 


Approved for public release; distribution is unlimited 
A NONLINEAR WAVE SHOALING MODEL 
FOR ALONGSHORE VARYING BATHYMETRY 
David M. Ruth 


Licutenant Commander, United States Navy 
B.5., United States Naval Academy, 1991 


Submitted in partial fulfillment of the 
requirements for the degree of 


MASTER OF SCIENCE IN METEOROLOGY AND PHYSICAL 
OCEANOGRAPHY 


from the 


NAVAL POSTGRADUATE SCHOOL 
September 2001 


Author: CAS /L 


David M. Ruth 


To 
Approved by: ae 


Thomas H.C. Herbers, Thesis Advisor 


Ethrand 8 far ofe 


Edward B. Thomton, Second Reader 


Mary L. Batteen, Chairman 
Department of Oceanography 


ill 


THIS PAGE INTENTIONALLY LEFT BLANK 


iv 


ABSTRACT 


This thesis proposes an improvement to present near-shore wave prediction 
models. Using weakly dispersive Boussinesq theory, the shoaling of directionally spread 
surface gravity waves over a beach with gentle gradients in the cross-shore and 
alongshore directions is examined. Following Herbers and Burton (1997), the governing 
fluid flow equations are expanded to third order and depth-integrated over the water 
column. A resulting amplitude evolution equation for a spectrum of waves is derived, 
which is the main result of this paper. New terms in the higher order result include 
effects due to alongshore bottom slope, higher order cross-shore depth variations, and 
non-linear quartet interactions. The linear terms in this equation are verified by analytical 
methods using linear finite depth theory. Example computations for a monochromatic 
wave train over a plane beach quantify some of the improvements of this result over the 
lower order model. Opportunities for further development and verification of this result 
are proposed, and recommendations for application of the result in its present form are 


outlined. 
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I. INTRODUCTION 


This thesis proposes an improvement to present near-shore wave prediction 
models. Accurate prediction of wave shoaling transformation is crucial for the successful 
planning of amphibious operations, providing a strong military motivation for this 
research. Additionally, the action of waves and wave-driven currents on a beach drives 


sediment transport and beach erosion, stimulating civilian interest in this field as well. 


As ocean surface waves shoal from deep to shallow water, amplitudes increase, 
wavelengths decrease, and propagation directions refract toward normal incidence to the 
beach as a result of linear processes. Additionally, non-linear interactions can cause 
significant energy transfers to wave components with both higher and lower frequencies 


than the original waveform. In deep water (Kh >>1, where « is the wave number and h 


is the water depth) and intermediate depths (Kh = O(1) ), nonlinear effects on waves are 


usually evaluated with finite depth theory on the basis of Stokes perturbation expansion 
for small wave slopes (ak <1, where a is the wave amplitude) (e.g., Phillips, 1960; 
Hasselmann, 1962). Nonlinear interaction between two primary wave components 
excites secondary waves with frequencies and wave numbers that are the sum (or 
difference) of those of the primary waves. In deep and intermediate water depths, these 
"triad" interactions are nonresonant, and the resulting secondary wave amplitudes are 
small compared to the primary wave amplitudes (Phillips, 1960; Herbers et al., 1992, 
1994). Similarly, tertiary waves may be excited by non-linear interaction between three 
primary wave components, constituting "quartet" interactions. For certain combinations 


of wave numbers quartet interactions are resonant, causing a gradual exchange of energy 
1 


between the wave components that is important to the evolution of wind-wave spectra 


over large distances (Phillips, 1960; Hasselmann, 1962; Hasselmann et al., 1973). 


Unlike the deep and intermediate water regimes, in shallow water (kh«1) the 
secondary response is strongly amplified due to near-resonance. The conventional 
approach to describe these interactions is the use of Boussinesq-type equations, which are 
based on the assumption that measures of nonlinearity, ¢ (=a/h), and dispersion 
a°(=(kh)’) are both small and of the same order. This approach reduces the three- 
dimensional fluid motion problem into a two-dimensional one by introducing a 
polynomial approximation of the vertical structure of the flow field, while accounting for 
weak non-hydrostatic effects due to the vertical fluid acceleration (see Madsen & 
Schaffer, 1998, for a comprehensive review). Boussinesq’s original work (1872) was 
restricted to a horizontal bottom; since then that work has been extended to an uneven 
bottom in one dimension (Mei & LeMéhauté, 1966), and later two dimensions 


(Peregrine, 1967). 


Herbers and Burton (1997, referred to as HB in the following) formulate shoaling 
evolution equations for the spectral amplitudes and phases of a directionally spread wave 
field propagating over a shallow, gently sloping beach with straight and parallel depth 
contours that are accurate to the lowest order (¢) in nonlinearity and dispersion (see also 
Freilich & Guza, 1984; and Elgar & Guza, 1985). Limitations of HB include: 1) neglect 
of alongshore depth variation, 2) restriction to gently sloping beaches, 3) restriction 
to small oblique wave angles, 4) a low order approximation of dispersion effects which 
limits the model to very shallow water (nominally depths less than 10m), and 5) neglect 


of nonlinear interactions between more than three wave components. In this paper, we 
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propose improvements that address the limitations mentioned above by carrying out the 
HB derivation to one higher order (<7) in nonlinearity, dispersion, and bottom slope, and 
allowing for weak alongshore depth variations. This higher order model still assumes a 
gently sloping beach and small oblique angles, but improves the accuracy of the model 
for finite bottom slopes and angles. Dispersion, directional spreading, and strength of 
triad interactions are all represented to one higher order, as well as the effects of the 
cross-shore bottom slope. New effects introduced at the lowest order are alongshore 


depth variation and nonlinear quartet interactions. 


An overview of the technique used to obtain this result is as follows: The first 
order horizontal flow and pressure is shown to be depth independent. The continuity 
equation is then depth-integrated to obtain a lowest order expression for the vertical flow. 
This vertical flow solution is substituted into the irrotationality condition and the vertical 
momentum equation, which are then depth-integrated to obtain the next higher order 
expressions for the horizontal flow and pressure, respectively. This process is repeated 
until third order expressions for the horizontal flow and pressure are obtained. Following 
Freilich and Guza (1984) and HB, an assumed first order spectral wave field is 
substituted in the governing equation for the depth averaged velocity potential, and the 
solvability condition yields a coupled set of first order equations for the evolution of the 


spectral amplitudes. 


A significant challenge to arriving at such higher order solutions is that the 
complexity of the polynomials at each level of the derivation grows significantly as 
higher orders of non-linearity and dispersion are considered. Madsen & Schaffer (1998), 
Agnon et al. (1999), and Gobbi et al. (2000) derive a variety of Boussinesq-type 
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equations to high orders of accuracy, from which this complexity is quite evident. 
Although equations derived in these earlier studies are more general than those that 
follow in this thesis, they are not cast in a spectral form and require prohibitively 


expensive full time-space domain solution techniques. 


This thesis is organized as follows. In Chapter II the formulation of governing 
fluid flow equations are expanded to third order and integrated over the water column. A 
coupled set of amplitude evolution equations for a spectrum of waves, the main result of 
this paper, is derived in Chapter III. In Chapter IV, linear terms in these equations are 
verified by analytical methods, and example computations for a plane beach quantify 
some of the improvements of this model over the lower order HB model. Finally, 


Chapter V summarizes these results, and identifies opportunities for further research. 


Hf. THIRD ORDER GOVERNING EQUATIONS 


A. ASSUMPTIONS 


The equations of motion and the continuity equation may be written in 


nondimensionalized form (neglecting viscous effects): 


Ou Ou 
—+(ueV)u+w—-+4 Vp =0, 
rr awa 
Ow Ow Op 
— +(uwV)w+ w—+—+1=0, l 
rr aki ars ) 
Ow 
Veu+—=0, 
: Oz 


where f is time, V denotes the two-dimensional gradient operator (O/0Ox,0/ Oy), u is the 
vector (u,v) of the horizontal (x,y) velocity components, w is the vertical (z) velocity 
component, and p is pressure. The x-axis points onshore, y points alongshore, and z 
points upward with z=0 corresponding to the mean sea surface. The surface and 


bottom are defined by 7(x, y,t) and z = —h(x, y), as shown in Figure 2.1. 


z= HELA 


Figure 2.1. Definitional sketch of variables and coordinate frame. Directionally spread 
waves propagate over a beach with gentle cross-shore and alongshore depth gradients. 
The variable @ denotes the wave propagation direction, k is the wavenumber vector, and 
h, is arepresentative water depth. 
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All variables in equation (Dare normalized using gravity g, density of seawater 


p, and a representative water depth /,. Irrotational flow is assumed, so that 


Ou _ Ow Ov Ow Ou Ov 


Sa Sat See. (2) 
Oz Ox Oz Oy Oy Ox 

The surface and bottom boundary conditions are given by 

On 

w=—+uVn, z=, 
at ) 1) 
p=0, z=", (3) 
w= —uwVh, Z= —h, 


The variables 7, u, v, w, and p are expanded in terms of the non-linearity 
parameter ¢, defined to be the ratio of wave amplitude, a, to water depth, h. In shallow 


water, v and w are both O(Kh) smaller than u (where « is the wavenumber), and thus the 


perturbation expansions of the scaled variables are given by 


n=en ten ten +. 
u = eu, + CU, + eu, sass 
veo(e, tev, tev; +...) (4) 
w=oa(ew, te°w, +e°w, +...), 


p=-ztepte’p,tept. 
where o = kh is a measure of dispersion. The independent variables are then scaled 


Din Wg? gi - (e % 


CSS 5 = = 75 9 
Ox Ox! Oy Oy! Oz Oz! Ot or 


so that the order of the terms appears explicitly when equations (4)]and (5)Jare substituted 
into the governing equations. In the classic Boussinesq approximation, dispersion o” 


and nonlinearity €¢ are assumed to be of the same order. The water depth is taken to be a 
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function of the slow variables (x, y) such that se gk and aE = 6 ky where 
Ox Ox Oy Oy 


B=WVh/kKh is assumed to be O(e). In this approximation, alongshore depth variations 


are O(«h) smaller than cross-shore variations. 


B. VERTICAL STRUCTURE OF THE FLOW AND PRESSURE FIELD 


In this section the irrotationality conditions, vertical momentum equation, and 


continuity equation are integrated over the water column at successive orders of € to 


determine the vertical structure of the flow and pressure field accurate to O(c’). First, 


we expand the irrotationality conditions (2) Jand vertical momentum equation (1b) using 


(4b-e) and Lsrlto O(e*): 




















Ou, Ou, Ow, | 2 ce Ow, | 3 
ak = aa —— a, = O ; 6 
Oz! es ev) lee eel (9) 
rr 4 [Oe OM | 62 OM _ OM | _ Ve), (7) 
Oz Oz = =Oy O Oy 
O O O O O O 
LS) eg | | ee? | eS Oe); (8) 
Oy Ox Oy Ox Oy Ox 
Op Op, . Ow Op; . Ow Ow. Ow. 
oP (22 + mh |, ~ ai +U, Fel +w, Ae! = O(e°). (9) 


Note that w), vi, and p; are independent of z’ to O(e). Substitution of the dynamic 


surface boundary condition (3b) in (4e) yields 


P; =": (10) 


Next we expand the continuity equation (1c) and the kinematic bottom boundary 


condition (3c) and collect terms of like order, 


Ou, , Ow, 
Ly 











Ox’ Oz’ =0, WwW, =0 at z'=—h, (11) 
Fr =O, ma ate 2 


At the lowest order, equation (11) integrates to 
w, =—(2/+h)—. (13) 


Substitution of w; and p; into and (9) and integration over the water column 


yields the vertical distribution of the second-order velocity and pressure field 











zZ Oru 

1 Uv te , (14) 
Zz” Oru 

v,=V, [Fre Ss. (15) 
Zz” Oru 

P2 =2-[E +i] 2, (16) 


where U>, V2, and P2 are integration constants. 


The constant P2 follows from substituting p; and p2 in and evaluating the 


dynamic surface boundary condition (3b) to O(e*) 


P,= tM i?) 


The remaining constants U2 and V2 can be expressed in terms of the second-order depth- 
averaged flow variables Uy and Vy (equations (23ylanalzelin Appendix A). With these 


results, the second order wave motion can be written as 








== 12 2) 
win [ vnc EO, (18) 
12 2 2 
—-. Zz / h O u, 
oe | 2 7 3 | A 
ae Oru 
Do — No = E+ n' Ayai ‘ (20) 


Having determined the second order pressure and horizontal flow variables, we 


may proceed as in the case of w; and integrate|(12)|to determine wo: 


2 (21) 


du, . Ov, 
ai 1 
Ox 


Ox’ — Oy’ 








ze hz’? — hh?z! | Ou 
n=[5 2 2 | a) 








Next we evaluate u3 by substituting wi, u2, v2, and w2 into equation|(6)|and integrating the 


Ole") terms with respect to z’: 











12 Dae 2 14 13 B12 4 
aid (ac dna 
x x Oy Xx (22) 
spe 
Ox Ox! 


Expressing the integration constant U3 in terms of the depth-averaged third order flow U; 


using equation (74)|(Appendix A), we obtain 





_ eS hz =z? — | Ou, 
























































se — 
BS aT 6 | 6 45S Ox" 
ge W\(Pu av 
—|—+hz' + 2 23 
c % " Ox'* — Ox'dy' ©) 
Oh Ou, . hn, Oru 
(Fe £ |) — 1 
CP aia aor 
Similarly for v3, 
12 2 2 14 13 212 4 
v; =V;— = + ha! ae ul = - ie = au 1 
2 Ox'Oy' Oy 24 6 6 Ox Oy 
(24) 
=e Oh Ou, . Oh Ou, 
Oy Ox' Ax dy" |’ 
and substituting for V3 using equation (Appendix A) yields 
yew eg 2 We? dS) tu (2? | Ou, | Oy 
enn 6 6 45}Ox"dy’ | 2 3 Jldx'dy' dy” (25) 
[244 Ou, Oh , Oy Oh), mh 9 
2)\ Oy’ Ox Ax’ Oy 3 Ox'dy' 








Finally, p3 is determined by substituting pi, p2, wi, and w2 into the vertical 


momentum equation (1b) and depth-integrating the O (c”) terms: 











a Vu, av au, (Ou,\ 
—P4+/—+h / 2 1 abel 
aa | 2 ‘| ax'or’ | dy'ar ax? lax! 


ghee re || Ou , Oh Ou, 
= + Bago ae 7° 
24 6 6 jOxOt Ox Ot 








(26) 





From equation (3b) it follows that 
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P\,_, =—E(m tem +e’n, + Ole*))+en 
2 Oru 2 3 4 (27) 
+e" |n, +enh—, + O(c’) | +2 (PR + O(e)) + Ole") = 0. 
Ox’ Ot 
Hence, 
Oru 
B=; mn Axford! (28) 
and 
z!? Ou, dv au, (du,\) 
= + —+h / 2 i Le 
Pas | 2 ‘| ox'at’ | By'at! || ax” ea 
(29) 








= zt het Wz”) dtu, oh Ou, | Ou, 
2° 6 6 \ax®ar | dxot ” axtat” 


In summary, the flow and pressure fields are described by equations (10) Jand 
13) Equations (10)}and k13)H20)Jare identical to those in HB. The higher order 
equations, improve the HB results by expressing the vertical flow structure as a 
fourth order polynomial in depth, and by including nonlinear and bottom slope terms that 
are neglected at lower orders. 


C. REDUCTION OF THE GOVERNING EQUATIONS 


Here the velocity and pressure distributions derived in the previous section are 
substituted into the governing equations to obtain depth-integrated horizontal momentum 
and continuity equations. First, the perturbation expansions (4) ]and scaling relations 
are substituted into the horizontal momentum equation (la). Collecting terms of like 


order, 
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Ou, . Op, Ou, Ou, . Op, | 
+ e|2 4 u, b+ 
ot’ Ox’ | ot’ a Ox’ — Ox" 
0 O(uu O 0 0 
lees te), Oi yg, Ste 4 BPs |= Ge 
and 
Ov, . Op, | OV, Ov, . Op, 
+P + 
bf Oy ak Ok On 





ies | Ov, Ov, Ov, ee Ov, | Op,|__ 
rs) / 


Substituting w;, u2, V2, P2, U3, v3 and p3 (equations (20) (23)}|(25)} |(29)) yields 


after some algebraic manipulation 








du, . On, Wl Bu, ay, 
-e pe es 
ot Ox’ — 3. | Ax” Ot’ ~— Ox'dy'dt' 


ho Oru O Ou 
3 oF fn Se] + Solaie) +» aa 


h? |. du, Ou Oru O Oru 
- 3 2 De! ax? wie a au 2 ? ote’), 


nh Ort, _ 7, Oh Oru, 
45 Ox"dt' — Ax Ox'dt' 


























(30) 











Ox'Ot 





dv, On, [dv, On, Ov, Wh Gu, 
pine | 7 rigs a eae Ta lait 
Ot Oy Ot Oy Ox’ 3 Ox'dy Ot 




















[Ov, On, Wl Fu, it Ov, | ht Ou {Oh Ou, oh au, 
ot’ dy’ —- 3. | Ax'dy'at' dy" at' | 45 Ox'dy'at' 2 dy Ox'dt' Ox Oy'dt" G1) 
ho Ou +u ov. 5 Mie, ov, 
20H Oeey| ae ae al 





Ou, O7u hr a Oru O Oru 
es ax'oy! 3 mae vn {n U||-o 


Finally, integrating the continuity equation (lc) over the entire water column, 
applying Leibniz Theorem, and the kinematic surface and bottom boundary conditions 


(3a,c), 


On oO a | 
ens ude + =O: 32 
Ot Oxd yi d Ye O) 


—h 
Substituting the expansions (4a-c) and|[(5)lyields 


aT 
nea tent Oe) 7) late Wie Ol) a 
; zi (33) 
2 O 2: 3 a 
+o ays vente v, + O(E ))dz = 0. 


0 
yi 


Now evaluate the integrals for each of the flow variables using equations (4a), (19) 


[53] ana BS] 


1 
f mae! = (ht+en, +n, )u, +O(e), (34) 


—h 











n 
i undz' =(h+en)u, —enh ot 1 4 O(e?), (35) 
—h 
f — nh? & 
i u,dz' = hu, + fe re =. 4O(e), (36) 
—h 
tT 
f v,dz! =(h+en,)v, + Ole”), (37) 
—h 
7] = 
Jf v,dz' = hv, + O(e). (38) 


—h 
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Combining and ordering the terms, 


On, Ou 
+h 
ot’ Ox Es 


| Om 
“Var! 


Ou, 
Ox’ — Oy" 


+h 














Ons 
Or’ 


du, Ov, 


+e 
Ox’ — Oy" 














+h (39) 














The system of equations and constitute the full set of governing 


equations to O (<’) . This system is now cross-differentiated to eliminate 7,, which gives 













































































Oru O ak on) 
2 
ar axl ox 
Puy O{,[du , an))__2(h (ru) 
Or = Ox Ox Oy Ox|3 Ox°dt? 
a (dh O( Ou,) 8 (mu) 
pat laser (Gaye (caceeaee aa 
~ Ox lax Ot\ Ox OX 
Qe: a oe 2 47 4 2 6 
Seg] Ca Oy | Rs OMe Lae cae Ov», | AL ui 
Ot Ox Ox Oy 0x00? OxOyat? 45 ax‘ar 
O|Oh— | Oh oO ne u.) ho Au, 
u, +—=v,|—ho— ae hi 
~ ax lax Oy OxOt ” xa SOF | "0x" 
Reca Ou, > Ou, aaa. du, ra a (Se) 
atl | dx * ax} atl | Ay J 2 Axdt\ Ax 
hr oO Oru Oo O 3 
Eh ma aye (hits + mat) Baa (= oe) (40) 





and 





Ov, O [i Ou, | 
















































































ar dy\. ax 
te avy, 8 Ou, Ov, at h O° (hu) 
or dy Ox Oy Oy|3 dxdt 
med eres eae 
Ov \ Ox Ot\ Ox OxOy 
2f@%_ 2,26, )|_ (ate, ay) ou 
or Oy Ox Oy 3 |dxdv0r Oy’ Ot? | 45 Ox’ Aydt? 
_d (ak an), Fu ha ay 
Oylax ? | dy || dydtl " oxat} 3 or" dxay 
9, % 7 or), O(, a ,# & eal 
atl | Ox 7° dx} Atl | Ay} 2 Oyat\ ax 














a ayael On| oxy ee 


where the primes are dropped from here on to simplify the notation. 
D. GOVERNING EQUATIONS IN VELOCITY POTENTIAL FORM 
Finally, expressing the horizontal flow variables as gradients of velocity potential 


functions @, (see Appendix B) allows equations and to be converted to the 


gradient form 








vf z9 h ad, 


















































Ot? Ox? 
ad, (06, &d,) W O*¢ 
+ ¢|—3 —-h| —3 + 3 |--—-— 2 
“lar ay] 3 dear 
~ 226 19 (2A 2, #4) 
Ox Ox 2 0t\ Ox Ox\" Ox 
ao, ad, a, h? fala ad, Od h* Od 
42/23 _4 cB mse’ 5] [ene ene 24 i eee 1 
* lar ve dy | 30° lar ay | 450xar 
_ Oh Od, _ Oh 0g, _, dh 0 1 Af A, 
dx Ox dy Oy ax dxor  Ot\" axar 
, 9 (89, 96,) ,18(A6) ,# a(aH) Wa (a4, ao, 
Ot| Ox Ox | 2 Ot Oy 2 Ot| Ax? 3 Ot\ Ox Ax’ 
d{ 6, O64) A{ a4, 3 
- + = +O = 0. 
Ox aC Ox mn Ox} dy 2 Oy (e ) (42) 





Since equation must hold in the entire horizontal plane, it follows that the 
expression inside the gradient operator (discarding a uniform time dependence) must 


vanish. Therefore, the final form of the governing equation is 


ad, = h ad, 
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se 





























Ox? 
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or Ox? Oy" 3 Ox’ ar 
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Ox Ox 2 0t\ Ox Ox Ox 
Oo, Fd, | Ob, _W a (io, o,)_h* 0%, 
Or Ox? Oy” 3 Of | Ox’ ~~ Ay’ | 45 Ox’ Or? 











dh 0b, Oh dG _, h_Oo, id A 


Ox Ox Oy Oy Ox Ox Atl" DAT 











d(d¢, 0¢,). 1 (4, es a (a4, ie d (dd, A°4, 
Ot| Ox Ox 2 Ot\ Oy 2 Ot| Ox? 3 Ot| Ox Ox? 
a( a0, 0¢,| O(_ d¢, 3 
= pac in peewee eels —0 
Ox G ax ax | yl Oy | (©) (43) 
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Hil. THIRD ORDER AMPLITUDE EVOLUTION EQUATION 


A. ASSUMED SPECTRAL WAVEFORM 


In the previous chapter, we derived a governing equation (43)] for the depth- 


averaged velocity potential function to Ole ). We now derive spectral amplitude 


evolution equations by assuming the lowest order wave field (7,¢,) to be a linear 


superposition of nearly plane, shoreward-propagating waves 


= S S54, (x,y) expliv, , (x,y,0)), (44) 
ies 
= ap Ee ci, (x,y,0)), (45) 


where w,=pAw and /,=qAl are the scaled frequency and initial alongshore 


wavenumber, with Aw and AI the separation of adjacent bands in the Fourier 


representation. The phase function y,, contains the fast variations of the lowest-order 


wave field, 
pg (% Dat) =\ a Tepe th eet (46) 


The complex amplitude A,, is a slow function of x and y 


OA, 
—£—=cR, +R 


1p.q 2p.4 


e . +0(e), (47) 


lp.q 


+0(e),  —t=es 


where Ripg, Rog, and Sj, account for amplitude and phase changes resulting from linear 


(e.g., dispersion, shoaling, refraction) and nonlinear (triad and quartet interactions) 
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processes, and changes in /, associated with alongshore depth variations are absorbed in 
the slow alongshore variations of 4,,. It will be shown that 04, /oy terms appear at 
sufficiently high order that they do not require explicit expression beyond O(e) : 


B. _ O(1) PROBLEM 


Define the linear operator £ as 


O° or 


It is easily verified that ¢, is a solution to 


£(4)=O(€). (49) 
To determine solutions to higher order problems, the residual right hand side of (49)|must 
be evaluated. To O(e" )s iL (¢,) is given by equation lonlin Appendix C. 
C. O(c) PROBLEM 


At O (e) , applying the results of Appendix C, 





BEE 





Oh —hw* A hl’ A 
2K hl +A,, naa PP» + q PY 





3 8) 
(50) 
3 


Ww 
+ a DL Ana Ayman 





expliy, ,| + O(e). 





Resonant growth of db, is prevented by the solubility constraint that the right-hand side 
forcing terms in |(50)|do not contain any free wave solutions of the general form 


F (x/ hi? 27): Since all these terms obey to O(e) the shallow water dispersion relation, 
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it follows that the right hand side of [eo must vanish. Solving for Rip, yields the lowest 


order term in the order amplitude evolution equation: 


Ring = Liy.gAp.g + 0, ys > Ay Ay 5 .y3 (5 1) 
where 


_ at | =| 


1d Ah Ox 6 2u, 
_ 3w 
Phen Eg 


This result was derived in HB. 
D. O(<”) PROBLEM 

To determine the remaining higher order amplitude evolution terms R», and Sip, 
we proceed in the same fashion as for Rip. Since L(¢,)=O(¢), , represents (to 
lowest order) free shallow water wave motion that can be absorbed in ¢,., and thus ¢, 
can be set equal to zero. Note, however, that 1, (evaluated in Appendix D) contains 
deviations from free shallow water wave motion that contribute to the forcing of ¢,. 
Tedious but straightforward evaluation of all the O (<”) forcing terms in equation [43)] 


(see Appendix C) yields 
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2P.4q 1p.q ee - 
Ww, Ox 3 WwW, Ox 


2hl ih OR 2hnvur ij 
Se APU. as +H ORes «|e Lot R 
P q Ww, 








ae (52) 


" Sh'’w, Oh 1, Oh 
45 2 


6 Ox w, Oy 








{=P ie a 










































































OD de a ace 
ip 2) 3), 2 = Ann Asn rmg-s-n feXPlivp,4|+O(e). 
Finally, substituting Rj, 4 (equation [(51)) and 
ORiyg _|—1 0h rake wl, Oh, | 
Ox 4h ax 16\0x 36 6 us IP 
enue {He 2) Beale. 
2p D)2 = nt A RA oO le) 








in equation |(52)| and setting the forcing terms on the right hand side equal to zero to 
remove all secular terms, the following relation for the higher order amplitude gradients 


Ropq and Sip is obtained: 


Y2q 
Ryn + Ww = Sip.4 = Ly y.qApg a ye Oe A (53) 
Pp r Ss 
F is Cs 5% » > As Hy, i; eee 2 


where 
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Z i ae |e l, Oh 
eS 
ax 2h'?w, dy 


2 
ey 


























alae De olf, 1 Ae 1 @&h 
360 | 12 Sw) 32hw,ldx}  8h"?w, ax?’ 
0, =zis on, | 3 (% 4), 3 (wi 2 
2 pqs 32h Ax 16h" 3 Ww, Shi? 3 w 
= Sus + 8w ws (3w, —w,)B —wLl, 
12h? 2h ux 
Swe OW, 
2p = 16h? 


E. EVOLUTION EQUATION 


Combining |(47)}|(51)| and|(53)} our final result in amplitude form is 


OA h'?1, OA 
Pay ; By = (Ligg ie Li, As, + ye Ge Sy Onecare) AeA ip 


7s oss De Ss Ai A, 7, eae + O (€* ) ‘ 





Ox Ww, 


(54) 


The real part of the linear term on the right hand side is the amplitude growth of the given 
mode due to shoaling (i.e., changes in the group speed); the imaginary part represents 
slow phase changes owing to dispersion, directional spreading, and bottom slope effects. 
The double and quadruple sums contain the slow amplitude and phase changes of a mode 
resulting from nonlinear interactions of all possible triads and quartets (respectively) in 


which the mode participates. 


The new terms (relative to HB: L2p4, Qrp¢7,s, and C2,) introduce improved 
representations of dispersion, directional spreading, triad interactions, and the effects of 
the cross-shore bottom slope (all to one higher order of accuracy), as well as the lowest 


order effects of alongshore depth variations and nonlinear quartet interactions not 
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represented in HB. The verification of many of these new terms, and a quantitative 


investigation of the improvements gained by them, is given in the next chapter. 
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IV. VERIFICATION OF LINEARIZED MODEL FOR AN 
ALONGSHORE UNIFORM BEACH 


A. VELOCITY PROFILES 


From linear finite depth theory, the horizontal velocity profile (in the x- 


direction) can be expressed in normalized form as 





u cosh («(z+h)) _ rh eosh(r(z +/)) 


— 1 (55) 
e = [cosh «(2 +h) dz sinh (xh) 


Expanding the hyperbolic functions for small «/ and using the geometric expansion for 


the quotient in equation (55)|yields 


(z+hy KS (z+h) 








K 6 
; Ph ay + O(Kh) | 
we (xh) (nh) aye 
1 eee ae 
Kh t+ + O(Kh) 




















taal 





Collecting and combining like terms, 


+0(kh)’. (56) 





u 512 h? a{z Age’ hz? ht 
==1+K |—+hz+—|+k ase = 

u 2 3 24 «6 6 45 
The z-polynomials at each order in equation |(56)| are in exact agreement with the 


polynomial coefficients at the corresponding orders for u2, v2, u3, and v3 (equations 


and \25)h. Figure 4.1 compares vertical structure of the first (u,), second 
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(u, + Eu), and third (u, + €u, +°us) order approximations of u with the exact linear 


finite depth solution Poe a wave with period 10s in 20m water depth. The first order 
approximation is the familiar depth-independent shallow water solution for the horizontal 
flow. The second order approximation better represents the depth dependence of the 


horizontal flow, while the third order approximation virtually replicates the exact finite 


depth solution. 
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Figure 4.1. Vertical profiles of normalized horizontal current . Note that the first order 
estimate is depth independent. Higher order estimates converge toward the vertical 
structure predicted by finite depth theory. 


B. DISPERSION 
In linearized form , the amplitude evolution equation for an alongshore 


uniform beach (dh/dy =0, 04/dy =0) is given by 
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Oh 
Ox 








Ww 2 
Ww 
Dp 
wp WE o>( te u,b : ec ae 


6 w, 360 12 8u} 32h?w lax} 8h!?w, ax’ 








(57) 








+i 





a 


The imaginary part of the right hand side represents slow phase changes as the 





wave component evolves in the cross-shore direction, that is, corrections to the cross- 


shore wavenumber, k. From equation it follows that k for a wave with radial 


frequency w and alongshore wavenumber / in uniform depth is given to O(«h)" by 


WwW hoy 11h??? v her o Bul” 7 p24 


k= 
jh’? 6 360 Qw 12 Sw 





(58) 


Now consider the exact & in finite depth. The alongshore and cross-shore 


wavenumbers satisfy the relationship 

+P =k’. (59) 
Using the finite depth dispersion relation, the cross-shore wavenumber can be expressed 
as 


4 
Ww 


= | — 60 
tanh? (Kh) 0) 


Since / is O(«h) , equation (eoylcan be written as the following Taylor series expansion: 


Ww? 7 yp? tanh («h) 7 1* tanh? (kh) 


= O(Kh)’. 61 
tanh (KA) ple 8uw° 2 (rh) (61) 





2a 


Now, expanding the hyperbolic functions on the right hand side of equation |(61)|in terms 
of Kh to O(«h)’ , and then expressing each «kh term as a function of w and / using the 


expansion of the finite depth dispersion relation, 


Ww =o GP 2) +0(kh)’, (62) 





we obtain 


Ww hi 11h we? h?P hw? yer 
=— + + > — = 
hv? 6 360 Qw 12 Bu" 











(63) 


in agreement with equation GS Figure 4.2 shows the first (O(1)), second (O(Kh) ), 


and third (O(Kh)") order approximations of the cross-shore wavenumber compared to 


the exact finite depth solution for a normally incident wave with a period of 10s. Figure 


4.3 shows the same comparisons for an incidence angle of 30 degrees. 


28 


020 degrees P= O01 ble 





Exact finiig capt 

isl order apps ceimation 
O 2 one approximation j 
a Tn archer ap proelmeienn 






1 


Wiayenunmbes 


0.088 } 


0.06 


Ci 


O.02 = 








to 20 ao 40 =a) oOo vo Bo oo foo 


Figure 4.2. Cross-shore wavenumber vs. depth for the first, second, and third order 
approximations are compared to the exact solution. Results are shown for a 10s period 
wave at normal incidence. 
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Figure 4.3. Cross-shore wavenumber vs. depth for the first, second, and third order 


approximations are compared to the exact solution. Results are shown for a 10s period 
wave at a 30 degree angle. 


Though not a significant improvement for depths less than 20m where Kh <1, the third 
order approximation does show a notable improvement over the second order for a wide 
range of intermediate (40-100m) depths. Note that the third order approximation 
eventually diverges from the exact result in sufficiently deep depths; as expected from a 
Taylor series expansion for small «h, this divergence grows rapidly with depth compared 


to the second order result. 


Figure 4.4 shows the same approximations of cross-shore wavenumber k 
compared to the exact finite depth solution for a 10s wave in a fixed depth of 20m, versus 


incidence angles. 
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Figure 4.4. Cross-shore wavenumber k vs. incidence angle 6,. The first, second, and 


third order approximations are compared to the exact linear solution. Results are shown 
for a 10s period wave in 20m depth. 


The first order approximation does not account for directional effects and has 
large errors. Both the second and third order approximations accurately account for 
directionality. The third order approximation is within a few percent of the exact 
solution even for angles up to 60 degrees. 

C. SHOALING 
The conservation of energy for a wave train with radial frequency w and 


alongshore wavenumber / may be expressed as 


O 
5, (Eee) = 0, 
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or 


aE _ 
Ox 











1% | P 222 (64) 


e. Ox rae Kk OX 


where E = |A) is proportional to the wave energy, and c,, is the cross-shore component 
of the group speed c,. For arbitrary finite depth, the (non-dimensional) linear dispersion 


relation is given by w* =« tanh(«h), yielding the following analytical expressions for 














bs Oe and oe, 
Ox Ox 
Ow wil Kh 
=e 65 
“8K oft ~ 
OG | oe Qwh? OK i Ww 1 2Kh Oh (66) 
dx | 2x? tanh? (2h)sinh(2Kh)} Ox sinh(2«h){ — tanh(2«h)J Ox’ 
2 
eae ~—-—-— a (67) 
Ox sinh’ (2h) +Kh a 
Substituting equations (67)|into equation (64)]yields 
OF | & Sinh (2h) = 2K (Kh) 
Ox Ie 5 1g ; 
2 [Sinn (2h) + ch tanh (2h) | sinh (2h) af ch 
2 2 
(68) 
al " : ‘ef? Oh 
tanh (25s (2c) + (1? =) [F sinh (20h) + a 
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In order to compare the exact, linear energy equation with the Boussinesq form 


(sa) we expand equation ég\lin a Taylor series for small Kh to O(«h) , which yields 


OE _ 


Ox 


Oh E. 


Ox (69) 


1, (nh) +a] 


2 2 
Kl 


2h 2h 2h 








Now, since E =|Af = AA (where "*" denotes complex conjugation), the 


amplitude evolution equation |(57)|yields the following energy equation (dropping the p,q 


indices) 


OE _ 404, gp OA _ 


Reeds (70) 
Ox Ox Ox 








By observation, it is immediately seen that equation (69) jand equation (70) jagree 


to first order. To compare the O(c) terms, note that uw =Kh(1+O(e)) and 


a (1 of O(e)) , and thus the O(c) terms in equation Lzoylare in agreement with 
equation (69)|to O(«?). 


The improved energy conservation properties of the present model can be 


observed in Figures 4.5 and 4.6, which compare the second and third order 


ae ; 1 Oc,. 
approximations of the energy flux coefficient ——-—* 


Coe OX 


to the exact linear finite depth 


solution. Again, results are shown for a 10s wave at normal incidence, Figure 4.5, and 


for a 30 degree incidence angle, Figure 4.6. 
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Figure 4.5. Energy flux coefficient —— co vs. depth for the second and third order 
bi, OX 
approximations are compared to the exact solution. Results are shown for a 10s period 


wave at normal incidence. 
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0 
Figure 4.6. Energy flux coefficient - =e vs. depth for the second and third order 


Cy OX 


approximations are compared to the exact solution. Results are shown for a 10s period 
wave at a 30 degree angle. 


Note that the second order approximation diverges from the exact solution even in 
relatively shallow water (<20m), and that this divergence is even more extreme for 
oblique waves. The third order approximation shows significant improvement over the 


second order in both the normal and oblique cases. 


To illustrate the improved shoaling characteristics of the third order model, the 
amplitude evolution of a wave shoaling on a plane beach was predicted numerically using 
a simple first-order finite difference approximation of the amplitude evolution equation. 


Figure 4.7 compares the evolution of |A/Ao| (where Ao is the complex amplitude of the 
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incident deep-water wave) predicted by the second (HB) and third order model with the 


exact finite depth solution. 
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Figure 4.7. Shoaling of a 10s period wave normally incident on a plane beach with a 1:50 
slope. Topography varies in the cross-shore direction only. Top panel compares the 


second and third order model results to finite depth theory. Bottom panel shows the 
depth profile. 


The significant improvement of the third order model over the second order HB 
model is evident in the predicted amplitudes in shallow water. For this relatively deep 
(25m) model initialization, the second order model shows the expected overshoaling, 
overpredicting the wave amplitude in 1m depth by 28%; whereas the third order model is 


within 1% of the exact finite depth theory prediction. 
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V. SUMMARY AND FUTURE RESEARCH 


In general, wave transformation on beaches is well-described by one-dimensional 
models. However, the directionality of waves is of crucial importance to a variety of 
nearshore processes including infragravity motions, alongshore currents, and sediment 
transport. Additionally, alongshore varying bathymetry complicates nearshore wave 
evolution. In this thesis, Boussinesq-type equations for weakly nonlinear, weakly 


dispersive waves are derived for a beach with weak alongshore depth variations. 


We extend the derivation of Herbers and Burton (1997) to one higher order in 
nonlinearity, dispersion, and bottom slope to derive a coupled set of amplitude evolution 
equations for a spectrum of waves. The main result, equation more accurately 
represents dispersion, directional spreading, and the strength of triad interactions than the 
HB model, and extends the effects of cross-shore bottom slope to one higher order. New 
effects introduced at the lowest order are alongshore depth variations and nonlinear 
quartet interactions. Linear terms were verified by analytical methods, and example 
computations in uniform depth and for a plane beach quantify some of the improvements 
of this model over the lower order HB model. In particular, it is demonstrated that this 
higher order model better represents the vertical structure of the flow field, notably 
improves the approximation of the cross-shore wave number both in intermediate depths 
and for larger directional spreading angles, and more accurately accounts for energy 
conservation. For a typical swell period and model initialization in relatively deep 
(continental shelf) depths, it is shown that the overprediction of shoaling is reduced from 
O(30%) for the HB model to less than 1% for the present higher order model. 


ai 


Further work is needed to verify other terms in equation No simple 
analytical results exist for beaches with two-dimensional depth variations. Therefore, the 
term that accounts for alongshore depth variation may have to be verified numerically 
using an appropriate refraction model. Higher order cross-shore bottom slope terms in 
equation involving the square and the curvature of the bottom slope are perhaps 
more difficult to verify. However, Chu & Mei (1970) derived similar terms for a Stokes 


wave in arbitrary depth. Their expressions are rather complicated in terms of hyperbolic 
functions. An expansion of the Chu and Mei result for small Kh to O(«h)" should 


render terms that correspond to the cross-shore bottom slope terms derived in this work. 


Finally, the coupling coefficients for triad and quartet interactions remain to be 
verified. HB show that in the limit of small wave amplitude and bottom slope the steady 
solutions to their evolution equations exactly match the second order bound waves 
predicted by finite depth theory in the shallow water limit. It should be possible 
(although more algebraically tedious) to extend this approach to one higher order and 


validate both corrections to the triad coupling coefficient and the quartet coupling 


coefficient in equation|(54) 
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APPENDIX A 


The average horizontal flow (u, v) over the water column is given by 


“= al udz' = Eu, +e7u, teu, +O(e*), 
v= ee | vOz' = o(ey, +e7v, +e°v, + O(*)). 


—h 


(71) 


The reciprocal of the water column height can be expressed by geometric expansion as 


1 ae! 1) 2 ™ No 3 
(a a ee 2 +O(e ). 


The integral term of u is expanded (note uw; is independent of z’ ) 


u) n 


uae! = | (eu, +2'u, teu, + O(e*) az! 


—h —h 


= (a ten ten, + O(")) EU, 











0 Em 0 
+e f mae! + [u,d2" +e [ udz' + O(¢*). 
zi 0 = 
Therefore, 
= l 0 ” 0 n l 0 
u=eu, te?|— | u,Oz’|te?|— | wdc’ +—4u +— | u,dz'|4+ Ole 
1 id 2 ed 2 h 2|z=0 id 3 ( ) 














and depth-averaged flow variables can be defined to correspond with each order of u: 
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u, = > f u,0z' (72) 





z=0 





aa uy u 
Wf us02! el ; 
= 


Substituting (14) and (22)|into (72)|for u2 and u3, respectively, and integrating yields 











= h? o 
=U, +> (73) 
and 
a ae h’ d'u, Wlu, 8, 
: > 45 Ax’* — 3. | Ax” Ox'dy! (74) 
4g Oh Ou, v hn, Oru, _ An Ou, 





Ox Ox’ = 3 Ox” )— 3 OX’? 


Similarly it can be shown that the depth-averaged alongshore flow variables are given by 





v=, (75) 
ie 2 92 
Vy 2 = au , (76) 
3 OxOy 


and 


Pu, . Ov, 

+ 
Ox'dy' dy’ 
hn, Ou, _ hn, Ou, 

3 Ox’? = 3 Ox'dy"” 
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(77) 
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APPENDIX B 


Consider the set of velocity potential functions | a 


OG LEE, iy 
Ox‘ Oy 


make necessary substitutions to the depth-averaged governing equations, these functions 
must be expressed in a depth-averaged form. For convenience, define d, to be the 


average value of ¢, from the bottom to the still sea surface 


oe! y 
b= J dz. (78) 


Since uw; and v; are independent of z, so is ¢,. Applying Leibniz’s Integration Theorem, 


0 


— 1d, 00, (&-%(-A)) an 
W,=— [Shae x +€é F re 


—h 





(79) 


To eliminate the explicit boundary value, ¢, (—h), from equation recall 


Bernoulli’s equation 


2 2 2 
pe ENE ip (80) 


Substituting the perturbation expansions (4b-e) and the definitions of ¢,, and collecting 


terms of like order, we obtain 


 ersaaar Ae (81) 





af a == = : (82) 


Hence, using equation|(78)| 
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so that 


ad, a he Ou, uy; 
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i= _ Ou 
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Integrating in time, setting the arbitrary integration constant to zero, 


Thus, equation (79) simplifies to 


_ h? Ou 

— ¢,(—h) = —-——. 
$,—4,(-#)=—E 
7 a Ob hm Oh 





b= = 
Ox 6 Ox Ox 


Similarly it can be shown that 


0 hm Oh 


> Oy 6 Ox Oy 





Applying the same procedure to the next higher order yields 


— Od, _ mh ou, 
> Ax 3 Ox’ 





+ O(e), 


poh), 
> Ay 3. Oxdy 
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APPENDIX C 


Substitution of the assumed waveforms for ¢, (equation (44) and n, (equation 
las} and the solution for 1, (2, = 0) derived in Appendix D (equation los) into the 


governing equation|(43)/yields the following expressions for each of the terms: 
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APPENDIX D 


To determine 7,, consider the surface boundary condition 


w=—+u—+v—, Z=>7,. (92) 


Substituting the perturbation expansions (4a-d) and scaling relations into equation 


(92)|and collecting terms of like order yields 


= O(e’). (93) 
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Now, substitute the assumed lowest order wavefield (44)]and \(45)|(dropping the primes 


to simplify the notation) to express each term involving 7,, u,,andv, in terms of the 


spectral amplitudes A, , 
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Finally, integrating equation|(94)] with respect to time and setting the arbitrary integration 


constant to zero, 
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